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Reliable determination of vortex parameters from measurements of the microwave 

complex resistivity. 
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We discuss and propose a complete data treatment, in close contact to typical microwave ex- 
perimental data, in order to derive vortex parameters, such as pinning constant and viscous drag 
coefficient (also referred to as "vortex viscosity"), in a way as model-independent as possible. We 
show that many of the accepted models for the complex resistivity can be described by a single, 
very general analytical expression. Using typical measurements of real and imaginary resistivity as 
a function of the applied field, we show that, even for single-frequency measurements, it is always 
possible to obtain (a) estimates of viscous drag coefficient and pinning constant with well-defined 
upper and lower bounds and (b) quantitative information about thermal creep. It turns out that 
neglecting thermal creep, in particular and counterintuitively at low temperatures, might result in 
a severe overestimation of the viscous drag coefficient. We also discuss the impact of thermal creep 
on the determination of the pinning constant. The present results might lead to a reconsideration 
of several estimates of the vortex parameters. 

PACS numbers: 74.25.Nf, 74.25.Qt, 74.72.Jt, 74.72.Bk 



I. INTRODUCTION 

The electromagnetic response of High-T c superconduc- 
tors (HTCS) in the mixed state is of great interest for 
fundamental physics as well as for technological applica- 
tions (e.g. pickup coils for magnetic resonance imaging, 1 
filters 2 -). In particular, controlling vortex pinning is es- 
sential to reduce power dissipation and signal noise in 
devices, whereas intrinsic macroscopic quantities such as 
the viscous drag coefficient are intimately related to the 
electronic states of vortex cores, and thus yield a great 
deal of information about the latter^d It is then desir- 
able to obtain accurate and reliable determinations of the 
vortex parameters from experimental data. As a matter 
of fact, however, in HTCS the estimates of the vortex 
parameters span orders of magnitude, even when similar 
compounds are measured^ The HTCS vortex dynamics 
is notoriously very complexj&i^ Thus, the determination 
of the vortex parameters from the experimental data is 
not straightforward (to say the least). Among the others, 
vortex-vortex interactions, which give rise to non-local 
response^ 1 - and disorder— play a fundamental role. As a 
consequence, the response of the vortex system can attain 
very complex frequency dependencies even in the linear 
regime of small currents ] 9 ' 10 : 1 WV 4 In those cases, it is 
rather difficult to reliably determine (or even define) the 
vortex parameters. 

A substantial simplification comes by increasing the 
stimulus frequency to sufficiently high values, e.g. to 
the microwave range. In this case the amplitude of 
the vortex oscillations becomes so small 15 that the sys- 
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tern can be treated in the local, single vortex limit. 
In this scenario a single-vortex, mean-field approach 
is usually profitably exploited to determine the vortex 
parameters^ - — such as the viscous drag coefficient, 
the depinning frequency, the pinning constant. Mean 
field models have been often extended with the intro- 
duction of field-dependent effective vortex parameters in 
order to account for inhomogeneities and finite elasticity 
of the vortex systemi 23 i 26 i 34 i 36 i 37 Mean field models can 
also be generalized rather easily to include the screening 
effects and the interactions between the microwave cur- 
rent and moving vortices i 20 ' 21 ' 25 As a matter of fact, most 
experiments concerning the mixed-state electrodynamic 
properties in the microwave range make use of mean-field 
models for the interpretation of the data and for the de- 
termination of the physics of the vortex matter. 

Due to the complexity of microwave measurements, 
high sensitivity is usually obtained by using resonat- 
ing techniques^ while wideband measurements (e.g. 
Corbino disk ^ 51 ' 52 ' 53 bolometry") are confined to small 
temperature regions. Thus, one has to deal with single- 
frequency measurements. Discussion of such measure- 
ments in terms of vortex parameters is only seemingly 
straightforward. In fact, it is quite delicate because the 
choice of the specific model can deeply influence the fig- 
ures which can be extracted, since a different dynamics 
(i.e., frequency dependence) of the data with respect to 
the model adopted can significantly vary the estimates of 
the vortex parameters (a discussion, in a different con- 
text, is given in Refs 13311551 ). 

An example is given in Refl56l. which presents mea- 
surements on Nb films: the flux flow resistivity as calcu- 
lated from the data displays very different field depen- 
dencies whether the pinning contribution is taken into 
account or not. Another example, specific to HTCS, is 



discussed in the review of Ref[5f there, it is shown that 
estimates of the viscous drag coefficient can vary up to 
two orders of magnitudes whether pinning effects are in- 
cluded or not. In addition to those remarks, one has to 
note that in HTCS thermal creep is unavoidable, at least 
above ~ T c /2: without a proper data treatment, the use 
of even the simplest vortex models would pose serious 
questions on the reliability of the estimates of the vortex 
parameters. 

From the above discussion, it is clear that the reliabil- 
ity of the determination of the vortex parameters from 
experimental data, and an evaluation of their model- 
dependence, is an important issue before those param- 
eters can be compared to microscopic theories. The fact 
that in most models the number of parameters exceeds 
the number of independently measured observables fur- 
ther complicates the problem. 

In this paper we intend to address the issue of a cor- 
rect and reliable determination of the vortex parameters 
based on typical microwave measurements. In so doing, 
we will show first that many mean-field models can be re- 
duced to a universal expression for the vortex resistivity. 
We will use this result to discuss the model-independence 
of the vortex parameters. Second, by exploiting physical 
and algebraical properties of the models for the vortex 
complex resistivity, we will show that many additional 
information, including creep, can be extracted with re- 
spect to standard analysis. Finally, we will apply the 
novel, extended analysis to some illustrative example. 

The outline of this paper is as follows. In Section [TT] 
we briefly recall the main mean field models for the vor- 
tex motion, with emphasis on their common features and 
differences, and we rewrite the resulting expressions for 
the complex resistivity in a single formulation. Simple 
numerical examples arc used to illustrate the lack of ro- 
bustness of the derivation of the vortex parameters with 
respect to different models. In Section InTl selected data 
of the microwave magneto-resistivity in the mixed state 
of some HTCS thin films, Tl 2 Ba 2 CaCu20 8+K (TBCCO) 
and YBa2Cu307_2; (YBCO), are discussed in detail in 
the light of the previous Section, to exemplify the pro- 
posed extended analysis. The results arising from the 
application of different models are discussed. In Section 
IIVI we summarize the results and prospect future work. 
The Appendices report most of the details of the calcu- 
lations. 



II. THE MIXED STATE MICROWAVE 
COMPLEX RESISTIVITY 

In this Section we consider the complex vortex resistiv- 
ity Pvm as calculated within mean field models. We focus 
on the local response of the fluxon system. In this case, 
the surface impedance Z can be derived when needed by 
properly taking into account the geometry of the sam- 
ple and of the experiment £1 We specifically discuss the 
common case where the microwave currents flow in the 



isotropic plane [(a, b) plane, in cuprates] of a uniaxial 
anisotropic superconductor, and the static magnetic in- 
duction field B is perpendicular to the plane of isotropy 
and to the alternate currents. The vortex motion resis- 
tivity p vm relates the alternate current density J to the 
electric field E = Bxv induced by the moving vortices. 
When the vortex displacement is sufficiently small to 
ensure the validity of the local limit, e.g. with increasing 
the driving frequency, the vortex velocity v is customarily 
determined by writing down the dynamic equation for the 
balance of forces (per unit length) exerted on a single 
vortc-i 5 - 58 - 59 - 60 - 61 



??v + anri x v + VU = J x n<& a 



thermal 



(1) 



where fi is the unit vector along the vortex and <I>o is the 
flux quantum. The overall current density J exerts the 
Lorentz force J x n<f>o. Thermal fluctuations give rise 
to the stochastic force Fthermah responsible for fluxon 
jumps between pinning sites. 

Power dissipation of moving vortices is represented by 
the drag force ?yv, where the viscous drag coefficient r\ is 
related to the relaxation processes of the quasiparticles. 3 
Thus, it is a quantity essentially connected to the micro- 
scopic electronic state. The viscous drag coefficient 77, 
which is also commonly referred to as vortex viscosity, 
is not to be confused with the vortex fluid viscosity dis- 
cussed in Refl62i and there denoted with the same symbol 
•q. The symbol r\ is here chosen following the common use. 

The perpendicular (Hall) force on a moving vortex, 
a#n x v, is described by the Hall coefficient an- The 
Hall angle is Oh = arctan(a/f/r;). 

The effects of pinning are represented by the force 
— VU where U is the spatial function describing the pin- 
ning potential. 

In the harmonic regime e 1 "', at sufficiently high an- 
gular frequency lo only very small oscillations, around 
equilibrium positions, are involved. Thus, one approxi- 
mates VU ~ k p v/(ko), where k p is the pinning constant 
(also indicated as the Labusch parameter a/,) and v/(iw) 
is the vortex displacement. This force is the result of sin- 
gle interactions between pins and vortices, of collective 
interaction of the ensemble of pins and the vortex mat- 
ter, and of the fluxon system elasticity itself. As such, it 
is clear that it contains much of the physics of the vortex 
matter. 

In the left-hand side of Eq. (Q} the balance of the vari- 
ous forces acting on a moving vortex is strongly frequency 
dependent. An important characteristic frequency is 
the so-called (de)pinning angular frequency ui p = k p /rj, 
which marks the crossover between elastic motion, domi- 
nant at lower frequencies, from purely dissipative motion, 
arising at higher frequencies. 

The role of vortex mass is a longstanding issue ) 63 i 64 
that seems far from an accepted solution. Accordingly 
to most estimates^ in the following we will neglect the 
vortex mass for the microwave frequency range we are 
interested in. 



Equation (p} is the starting point for many models. 
Different expressions of p vm have been derived depend- 
ing on the specific physics incorporated in the various 
terms of Eq.JT]). As important examples, the inclusion 
of finite vortex elasticity and/or collective pinning phe- 
nomena through properly defined field dependent pinning 
parameter o 9 i 26 ' 36 i 37 yielded to a much wider application 
ofEq.p. 

In the following, we briefly recall several specific 
models for the vortex resistivity, and we show that all 
can be reduced to a single analytic expression. 

Gittleman and Rosenblum (GR) model. — In this sem- 
inal work^ no thermal and Hall terms were considered: 
Fthermai =0, a H =0 in Eq.©. Thus: 



Pvm,GR 



$> B 



V 



1— i=» 



(2) 



In this model r\ and u p can be directly calculated from 
the data by simple inversion. In the high frequency limit 
(u > Up) Pvm.GB. -* Pff, being p ff = $ B/r] the free 
flux flow resistivity. Equation (J2J) gave for many years the 
theoretical grounds for the interpretation of microwave 
and radiofrequency (rf) data, and it has served as an 
essential interpretative tool for many experiments per- 
formed in HTCS.i£i22^2&2LS&£LS£ 

We note here the relevance of the so-called r- 



parameter: 



_ 9(P„,Tl,Gfl) 



In general, it gives a measure 



' K(p„m,Gil)" 

of the relative weight of the reactive to resistive response, 
given essentially by the elastic and dissipative response 
of the vortex. However, in this specific model it directly 
yields r — u p /u. The peculiar, quantitative role that 
r plays in the GR model, together with the fact of 
being an experimental quantity, make it an important 
parameter in the discussion of the data (see also the 
thorough discussion in Rcf,69|). 

Coffey- Clem (CC) model. — By considering nonzero 
Fthermai, and assuming a sinusoidal pinning potential 
U(x), Coffey and Clem calculated the vortex resistivity 



as 



20 



Pvm,CC 



<&o^ + if- 



V 



1 



where e is a dimcnsionless creep factor. In the assump- 
tion of a uniform periodic pinning potential of height 
Uq and by defining a normalized energy barrier height 
u = U (T,B)/(K B T) (K B is the Boltzmann constant), 
the following expressions hold: 



Mu/2)]- 



to Eq.@, < e < 1. In the limit e — > 0, u — v u p and 
Eq.Q reverts to the simpler GR model. For e — » 1 
(Uq — > 0), as it can happen near the critical temperature 
T c , thermal creep completely washes out pinning so 
that p vm reverts to pure flux flow pff = ^qB/tj at all 
frequencies. 

Brandt (B) model. — Brandt considered the creep ef- 
fects by introducing a phenomenological thermally relax- 
ing pinning constant k p , t (t) — k p e~ t ^ Tr fRefl70l) so that 



the corresponding vortex resistivity is 
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Pvm,B 



1] 1 + 1LUT 



(6) 



where e' = 



T v +T r 



— takes the role of a creep parameter, 
: T p e u °/ KBT , being t p = l/u p the usual 

(de)pinning characteristic time and Uq the pinning 
potential barrier height. It is worth noting that, within 
this model, a tighter algebraical condition < e' < 1/2 
holds. This model does not require any specific as- 
sumption about the pinning potential landscape, thus 
allowing in principle to include also other mechanisms 
like quantum flux creep^ On the other hand, it cannot 
be used for too high creep rates (Uq — » <=> e' — > 0.5), 
since in this limit it does not recover the expected pff 
limit for any finite or zero frequency. 

Two-mode (TM) model — The GR, CC and B mod- 
els considered a single oscillatory mode for the vortex 
motion. A second dynamic mode of the fluxon system 
arises from the bending of the flux linesj 21 ' 72 This addi- 
tional mode strongly affects p vm in the case of surface 
pinnng, that adds up to the bulk pinning. In the present 
case of thin films, in absence of thermal activation the 
following expression takes placed 



Pv 



*oB 



1 



.,TM 



V i 



(7) 



where u Pte ff t TM = k p /r] + 2ei<&o/(lr)d) — ui p + ui ps is an 
overall pinning frequency combining the effects of bulk 
pinning (u p ) and surface pinning (u ps ), d is the thickness 
of the thin film, Ei is a (field dependent) fluxon line 

(3) energy,— and I is a characteristic length describing the 
surface pinning. 14 When surface pinning is negligible, 
I — ► oo. Equation ([7]) coincides with the GR expression, 
Eq.Q, once the pinning constant is redefined with the 
additional contribution given by surface pinning. By 
contrast, for bulk geometries (beyond the scope of this 
work) one has to take into account the bending of flux 
lines and the redistribution of the currents, so that the 

(4) full expression becomes rather complex<2L 



wo 



Vp 



1 h(u/2) 
l-eI (u/2) 



(5) 



where Iq and I\ are the modified Bessel functions of the 
first kind and orders and 1, respectively. According 



Universal expression. - - Despite their different ap- 
proaches, all the previously recalled models for p vm can 
be cast in the following single analytical expression: 



Pvva — Pvm.l \ iPvm,2 



<& Q B e e f f + iuT e ff 

r) ef f 1 + IWTeff 



(8) 



where r e /y is the main time constant governing the vortex 
oscillations (usually linked to the pinning characteristic 
frequency), and the dimensionless parameter e e // is a 
measure of the weight of creep phenomena. It is bound 
to the interval [0,1] in order to recover the correct zero 
frequency limit, which must be < p vm < pff. We note 
that, according to Eq.flSJ), p vm ,2 > always. 

Inclusion of the Hall term. — For completeness, we 
mention that also the Hall term can be incorporated in 
the models for p vm . Neglecting creep for the sake of 
simplicity, the straightforward extension of the GR model 
yields) 7 - 3 - 



Pvm,GRH 



$0# 



tan 



(9) 



-*H 



valid at constant microwave currents (it is useful to re- 
call that most of cavity/resonator experiments meet this 
condition, while e.g. Corbino disk experiments do not). 
The Hall term introduces a different frequency depen- 
dence: a Bode analysis reveals a pole in zero, as in the 
previously discussed models, and two additional complex 
poles which introduce a resonant response. Nevertheless, 
one can formally cast Eq. © in the form of the universal 
expression, Eq.®, with e e // = and defining an effec- 
tive viscous drag coefficient ?/ e //,_ff — V + a 2 /i] - often 
considered in microwave experiments^^ - and 



1 



1+ 21. 



r eff 



Wp,eff,H = Uj 



tan" Oh 



P (l + (^) 2 )(1 + 



taii. z e H ) 



(10) 



Wp,eff,H formally plays the role of an effective pinning 
frequency, but retains a dependence from the measur- 
ing frequency ui. For not too large Hall angles (i.e., for 
tan 8h < 7r/4, which is a condition generally expected 
in cuprate o 40 ! 43 ' 74 ), oj p , e ff y H is always a positive quan- 
tity. When interpreting single- frequency measurements, 
one has to bear in mind that the presence or absence of 
the Hall term might not be an obvious issue. At very 
high measuring frequency (in the THz range) one should 
take the Hall term into account^ 

The reduction of many models to the single analyti- 
cal expression, Eq.©, allows to conceptually divide the 
whole data analysis and interpretation process into two 
steps. 

First, using the sole assumption of standard vortex dy- 
namics, one can derive from the experimental data sev- 
eral model-independent quantities: T] e ff,r e ff,e e ff- Al- 
though their physical meaning can be fully determined 
only with the choice of a particular model, their deter- 
mination relies only upon Eq.©, and thus it applies in- 
differently to all the specific models described by Eq.©. 

Second, the full physical interpretation of the effective 
parameters is given after the choice of a specific vortex 
model. This last step can lead to quite different results, 
hence its delicate nature. 

Accuracy issues, applicability limits, general con- 
straints are thus interesting for the analysis of the 



experiments. 

Numerical examples. - - We now illustrate with two 
numerical examples the effect of neglecting flux creep in 
the interpretation of the data, and then the need for a 
more accurate data treatment. To do so, we first generate 
data for p vm (B) using complete models (e.g., the CC or 
B models) and reasonable (according to the literature) 
values of parameters. Then, we use a simplified model 
(GR) to evaluate ij and k p from the as-generated data. 
We show that this procedure, which mimics a very widely 
used approach to the interpretation of the experiments, 
can (1) yield wrong estimates for the vortex parameters 
and (2) drive toward a complete misinterpretation of the 
physics at the origin of experimental data. 

In the first example we compute p vm (B) — p vm ,i{B) + 
ipvm,2(B) at 25 GHz by means of the B model, selecting 
parameters appropriate for a large elastic response (as, 
e.g., in YBCO with strong pinning below 80 K): r\ = 10 
Ns/m 2 , k p = 3 x 10 4 N/m 2 and a small but finite creep 
factor, e' = 0.15. 




FIG. 1: Simulated p V m(B) vs field through the B model (n = 
10" 7 Ns/m 2 , k p = 3 x 10 4 N/m 2 , e' = 0.15) representing an 
YBCO with strong pinning, T ~80 K and «/=25 GHz. 

As seen in Fig[TJ both p vm ,i and p vm .2 are linearly 
proportional to the field B, which is usually taken as 
an indication for the absence of creep effects. However, 
applying the GR model to these same p vm (B) data one 
obtains r\ = 1.57 x 10~ 7 Ns/m 2 , largely different from the 
"true" value. By contrast, one would have k p = 2.7 x 10 4 
N/m 2 , 10% far from the true value. It is noteworthy that 
neglecting creep has the largest influence on the estimate 
of rj. This point will be further investigated later. 

In the second example we compute p vm (B) at 13.03 
GHz and 83.5 K by means of the CC model, using pa- 
rameters that were found to describe excellently real data 
taken on a YBCO thin film: 9 k p = 4.5 x 10 3 N/m 2 , 
t] = 3.5 x 10 -8 Ns/m 2 , and a field-dependent creep pa- 
rameter (upper panel of Fig(2j) calculated with a normal- 
ized pinning barrier u(B)/2 = 3.5 T/B. Applying the 
GR model (which neglects flux creep) to the as-generated 
data, one would find fictitious field dependences in 77 and 
k p (center and lower panels of Figure[2] open circles). We 
note that the fictitious field dependence of k p (B) is par- 
ticularly dangerous for a proper interpretation: although 
completely artificial, it is potentially verisimilar, since it 




the vortex motion contribution, p vm . Typical measure- 
ments of p vm {H) at similar reduced temperatures for 
both samples are shown in Figures [3] and [U The quan- 



FIG. 2: Vortex parameters used to reproduce real data of an 
YBCO thin film at 83.5 K and i/=13.03 GHz. Full symbols: 
values proposed in Ref 9 through CC model fits (k p = 4.5 x 10 3 
N/m 2 , 77 = 3.5 x 10~ 8 Ns/m 2 , u(B)/2 = 3.5 T/B). Open 
circles: values calculated through GR model. 



could be referred to a change in the pinning of the fluxon 
system i 59 i 75 Thus, the choice of the model has a very 
strong impact on the evaluation and interpretation of the 
vortex parameters. In the following Section we illustrate 
a method for the analysis of the data which alleviates this 
issue, giving reliable and possibly model-independent es- 
timates of the vortex parameters and, most important, 
estimates of the error bars. To keep contact with exper- 
iments, we illustrate the procedure making use of mea- 
surements of the complex magneto-resistivity on various 
HTCS samples. 



III. EXTENDED ANALYSIS OF 
EXPERIMENTAL DATA 
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FIG. 3: Vortex resistivity p vm vs applied field H in the 
TBCCO sample at T = 81.4 K (T/T c ~ 0.782). Upper panel: 
parameter r — p V m,a/pvm,i- Lower panel: p vm ,\ and p vm ,2- 

tity r — p vm ,2/pvm,i is also shown, given its relevance in 
the analysis. Due to unavoidable numerical uncertain- 
ties in the calculation of r at low fields, we will limit the 
analysis of the data at fields poH > 0.1 T. 

In both samples a significant magnitude of the reac- 
tive component can be observed: r > 1. In TBCCO, the 
upward curvature in the data (more evident in p vm ,\) 
suggests field dependent vortex parameters, with a pos- 
sibly relevant flux-creep. In YBCO p V m,i(H) is almost 
perfectly linear in H , consistent with a simpler scenario 
of constant viscous drag coefficient and negligible flux- 
creep. We note that p vm in TBCCO is much larger than 
in YBCO at the same fields, by about one order of mag- 
nitude. 



A. Experimental data 

The measurements here used for illustration were taken 
on a TBCCO thin film (T c ~ 104 K) and on a YBCO 
thin film (T c ~ 90 K) with BaZr03 inclusions. Details 
of the film preparation have been given elsewhere i 76 ! 77 
Microwave measurements were performed by means of a 
sapphire dielectric resonator— at v — 47.7 GHz, with a 
static magnetic field p,<yH ~ B perpendicular to the sur- 
face of the samples (aligned with the c axis) . The vari- 
ation of the microwave resistivity with the field yielded 



B. Conventional analysis 

We first apply the conventional approach (GR model, 
no flux-creep) to our data. By inverting Eq.([2|), tjgr 
and fc Pl G-R are directly obtained from the data (the sub- 
script "GR" stands for zero-creep derived values), and are 
reported in Figures and [6] for TBCCO and YBCO, re- 
spectively. In YBCO k Pt GR is constant and tjgr shows a 
weak field dependence (approximately within 15% of the 
average value). Absolute values are in agreement with 
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FIG. 4: Vortex resistivity p vm vs applied field -ff in the YBCO 
sample at T = 70.9 K (T/T c ~ 0.787). Upper panel: param- 
eter r = pvm,,2/pvm,l- Lower panel: p vm ,i and p V m,2- 
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FIG. 6: Vortex parameters as derived within the GR model 
in YBCO at 70.9 K. Main panel: viscous drag coefficient tjgr- 
Inset: pinning constant k Pi GR- 



shown, those features could well be an artifact of the sim- 
plified analysis. Since both features might represent non- 
standard vortex dynamics, we now illustrate a method to 
assess the reliability of the calculations of the vortex pa- 
rameters and, as a consequence, of the physics that can 
be derived from them. 
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FIG. 5: Vortex parameters as derived within the GR model 
in TBCCO at 81.4 K. Main panel: viscous drag coefficient 
rjcR- Inset: pinning constant k Pi GR- 



the literature (see Ref|5| and references therein). Apart 
from minor considerations, these results seem in line with 
conventional vortex dynamics, and would not stimulate 
particular comments. By contrast, in TBCCO ttjqr is ap- 
proximately field-independent but attains very small val- 
ues, about one order of magnitude smaller than typical 
values in YBCO. In addition, k p ,GR decreases markedly 
(by a factor ~3) with the field. However, as previously 



C. Extended analysis 

We first exploit the analytical and physical features of 
the universal expression, Eq.®. Only as a second step 
we make use of specific models. 

Creep factor — The creep factor e e ff cannot be di- 
rectly obtained from the data. However, it is possible to 
give an upper limit e e ff >ma x according to the following 
expression^ derived in Appendix A (we recall that r is 
an experimental quantity): 



e e // < e e ff,max(r) = 1 + 2r 2 - 2ry/l + r 2 



(11) 



This is an important result, that will be used throughout 
the remaining of this paper. 

As an example, e e ff,max(B) at T=81.4 K for the 
TBCCO sample is plotted in Figj7] (upper panel). It 
can be seen that e e ff iTnax is an increasing function of B. 
In addition, we have found that it is also an increasing 
function of T, consistently with the thermal origin of 
this parameter. 

Viscous drag coefficient — An important constraint on 
the viscous drag coefficient 77 -which we recall is also re- 
ferred to as vortex viscosity- can be derived (calculations 
are reported in Appendix A): r\ always falls in the range 



Vgr 



1 + £ e//, ; 



< V < VGR 



(12) 
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FIG. 7: Model dependence of vortex parameters in TBCCO 
at 81.4 K. Upper panel: maximum creep factor e e ff,max(H). 
Inset: corresponding minimum barrier energies u m i„ of the 
pinning potential within the CC (full dots) and B (open sym- 
bols) models. Lower panel: k p within the CC (maximum 
values, full dots), B (maximum values, open symbols) and 
GR (down-triangles) models. 



Note that Eq. (fT2")) gives model-independent (in the mean- 
ing of "all the dynamics that can be described by Eq.®" ) 
ranges in terms of experimental quantities: tjgr and 
teff.max are directly calculated from experimental data. 
The values obtained through the use of the GR model 
always represent an upper limit for 77. Eq. (|12[) gives a 
counterintuitive result: the GR estimate, tjcr, turns out 
to be a more precise approximation to 77 when e e ff tTna x 
is large, most likely at high temperature. In other words, 
in presence of a large reactive component -i.e. large r- 
the application of the GR model to low temperature data, 
where creep is thought to be small or neglibible, yields 
more uncertain estimates of n. In this case it is true 
that e e ff iTnax is small (and that the real e e // can be even 
smaller), but neglecting it altogether brings an error up 
to a factor of 2 on n. It is not excluded that the present 
result might bring a reconsideration of many of the values 
obtained for r\ in the literature. 



Thus, estimates for r\ can be accompanied by signifi- 
cant error bars. Errors can be reduced by noting that 
the uncertainty on the creep factor, < e e ff < e e ff. ma x, 
does not uniformly reflect on 77. In Appendix B we discuss 
in detail a statistical approach. The main result, that we 
state here, is that even in the worst case (77) ± 20% repre- 
sents a 90% confidence interval, where (77) is the expected 
value of 77. 

We stress again that the results derived up to now 
stem from the general expression, Eq.®, and thus are 
independent from the specific model adopted. In par- 
ticular, the experimental data here reported in TBCCO 
actually indicate a very small viscous drag coefficient^ 

Activation energy. — From the maximum creep factor 
eeff.max it is possible, in principle, to derive a minimum 
barrier energy t/o, m m of the pinning potential. To do so, 
a specific model must be chosen. The widely used CC and 
Brandt models usually yield similar results, in the range 
where the Brandt model can be applied {e e ff, max < 0.5, 
that is r > 0.35), as reported in the inset of the upper 
panel of FigJT] The main difference is a roughly vertical 
translation, arising from the different treatment of the 
thermal creep contribution. 

As an illustration, in TBCCO we obtain (7o, m j n (S=0.2 
T) ~20 meV and 40 meV in the CC and Brandt model, 
respectively. Similar values have been reported in 
YBCO, as determined through microwave multifrc- 
quency measurements ] 10 i 13 

Pinning constant. — The uncertainty in the determi- 
nation of k p behaves in the opposite way with respect to 
the viscous drag coefficient r\: the GR value k p ^R repre- 
sents the lower bound, whereas finite values of e e // yield 
larger k p , up to a maximum value k p-max which is model 
dependent. An interpretation of the data through the 
Brandt model yields k Pimax ~ k p ,GR, thus giving very 
small uncertainties (see comments in Appendix B). How- 
ever, an interpretation of the same data within the CC 
model yields a larger uncertainty, as represented in lower 
panel of Fig(7] Care must be taken before driving conclu- 
sions from an analysis of k p : the model chosen -differently 
from 77- plays a very important role in the derivation of 
k p from the data. Only the lower bound is model in- 
dependent. The uncertainty in k p has a trend opposite 
with respect to 77. As exemplified in FigO the relative 
uncertainty grows as k p decreases, which coincides with 
decreasing r. 

Coming back to the illustrative data reported here in 
TBCCO we observe that, whatever the model chosen 
(CC or B), the uncertainty band of allowed values is 
narrow enough to unambigously determine that k p is a 
decreasing function of the field. The study of the physics 
involved in this field dependence is beyond the scope of 
this work, and it will be investigated in the future. We 
only mention that such a behavior may originate from 
rather different physics. For instance, a finite, small 
density of effective pinning centers can yield a decreasing 



average pinning strength over isolated vortices/vortex- 
bundlcs as the increasing number of vortices with 
the field progressively exceeds the number of pins. 59 
Another possible scenario could be related to some 
field-driven transformation of the fluxon lattice, impact- 
ing vortex elasticity properties and, ultimately, pinning. 8 

Reanalysis of YBCO data — For completeness, we now 
consider the data^ 8 - taken on the YBCO sample reported 
in Fig[5]and analysed according to the conventional anal- 
ysis as reported in Fig[6j The extended analysis con- 
firms the results for the viscous drag coefficient: its be- 
haviour is fully compatible with a field independent value, 
and numerical values are in agreement with literature. 
The pinning-related parameters extracted through the 
extended analysis and reported in Fig[S] add some piece 
of information. The creep factor e e ff^ max and the nor- 
malized minimum barrier energy u m i n exhibit a similar, 
but weaker, field dependence as in TBCCO. Irrespec- 
tively of the choice of the CC or B model, k p is compat- 
ible with a field-independent behavior, in clear contrast 
with TBCCO. This fact suggests (a) the existence of pin- 
ning wells steeper in YBCO than in TBCCO and (b) the 
presence of additional field-induced effects involving the 
fluxon system in TBCCO. 



IV. SUMMARY 
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We have reconsidered the problem of the determination 
of the vortex parameters from measurements of the com- 
plex resistivity. We have reviewed several widely used 
mean- field models for the vortex dynamics, and we have 
shown that they can be cast in a single, general expres- 
sion. We have shown that, independently on the model, 
neglecting creep effects gives rise to a large uncertainty 
on the viscous drag coefficient. This property holds true 
in particular when creep is small (but nonzero). This 
finding might stimulate a reconsideration of the values 
for the viscous drag coefficient reported in literature. We 
have proposed an extended analysis which allows to eval- 
uate and keep under control the uncertainty inherent in 
the evaluation of the vortex parameters. We have also 
exploited some model-dependent implications. Overall, 
with respect to the conventional approaches the analy- 
sis here presented allows for the extraction of more, and 
more controllable, information from the data for the com- 
plex resistivity in the vortex state. 
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FIG. 8: Model dependence of vortex parameters in YBCO 
at 70.9 K. Upper panel: maximum creep factor e e ff,max(H). 
Inset: corresponding minimum barrier energies u m in of the 
pinning potential within the CC (full dots) and B (open sym- 
bols) models. Lower panel: k p within the CC (maximum 
values, full dots), B (maximum values, open symbols) and 
GR (down-triangles) models. 



V. APPENDIX A: CONSTRAINTS ON VORTEX 
PARAMETERS. 

In this Appendix we consider the universal Eq.([5]) for 
the vortex resistivity p vm '- using only very general ana- 
lytical and physical conditions, we derive relations and 
constraints among the physical quantities involved in the 
model. 

As a first step, we consider the following ratio of ex- 
perimental quantities in the light of Eq. (J5J) : 



Pvm, 2 
Pvm,l 



1-e 



eff 



1 



(13) 

W e e ff 

where w — l/(wr e //)- The ratio r belongs to the interval 
[0, oc]. By rearranging the above relation, one finds: 

{re eff )w 2 + (e e// 



1W + r = 



(14) 



The physical meaning requires vj > and 1 > e e ff > 
(see Section |TTJ) . Thus, from Eq. (fT4"|) one readily obtains 



l-€ cff -VA 
2re e ff 



(l-e e// -v / A)(l- £e// ) x I 7 ° ( ^ /2 ) i CC 

'-"■" I 1, B 



K (Fig(9j) shows that the expectation value (77) is closer 
to the zero-creep value, r]GRi than to the middle of the 
allowed interval, indicating that large e e // weakly con- 
tribute to (rj). In fact, plotting the 90% fidelity band 
(shaded area in Figj9j) starting from the zero creep value 
ijgr, one notices that the uncertainty on e e ff affects very 
unevenly the uncertainty of r\. Interestingly, this result 
relies only on the general expression, Eq . ([5]) . 
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TABLE I: Explicit expressions for the vortex parameters ac- 
cording to Equation (Jsj) as a function of e e ff and of the exper- 
imental quantities, p vm ,i and r. For the sake of compactness, 
we denned A = (1 — e e //) 2 — 4r 2 e e //. The pinning constant 
is easily calculated from the definition k p — ui p /rj. Here, Iq 
stands for the inverse function of Iq. 



the constraint: 



e e ff < £eff,max{r) = 1 + 2r 2 - 2ryl 



(15) 



which gives an upper bound for the creep factor that 
can be directly calculated from the data. The function 
£eff,max{r) decreases monotonously from 1 to when r 
spans the interval [0, oo]. The vortex parameters rj and w 
are correspondingly delimited within allowed ranges that 
can be easily determined through straightforward alge- 
bra using the explicit expressions for the various vortex 
parameters reported in Table |TJ 



VI. APPENDIX B: STATISTICAL ANALYSIS 

In Appendix A we derived the maximum variabil- 
ity ranges for various vortex parameters. However, 
more information can be gained by a statistical ap- 
proach. We ascribe the maximum uncertainty to e e //' 
for each (B,T) point e e // can be treated as a ran- 
dom variable with rectangular probability density in 
the [0, e e ff imax (B, T)] range. Thus, other quantities 
can be described in terms of derived probability den- 
sities, computable through standard theorems for the 
functions of random variables^ The expectation value 
for a generic quantity a derived from e e // is (a) = 

J e//,max a(e e ff)de e ff. The illustrative example 

of the viscous drag coefficient rj in TBCCO at T=81.4 



510 




^ H(T) 



FIG. 9: Estimates of rj as obtained from the same set of data 
in TBCCO at T=81.4 K. Full symbols: tjgr, assumption of 
zero creep. Open circles: expected value (77), assuming that 
the creep factor has flat probability in its allowed range. Con- 
tinuous line: minimum allowed value of rj, in the assumption 
that the creep factor attains the maximum value compatible 
with the data. Shaded area: 90% fidelity band, including the 
zero-creep value. 

Similar considerations can be done for the pinning con- 
stant k p , but in this case a specific model has to be cho- 
sen. We discuss some relevant cases. Within the Brandt 
model the uncertainty on k p = tOp/rj is vanishingly small. 
Numerically, this is a consequence of the almost exact 
cancellation of the uncertainties on 77 and u> p when r is 
sufficiently high. This fortunate combination does not 
take place in the CC model. To illustrate further this ef- 
fect, Fig llOl reports the maximum deviation k Vimax jk Vi GR 
as a function of r for both models. It is apparent that 
for r > 1 the Brandt and GR values are almost coinci- 
dent (kp^max/kp^GR — 1), while the uncertainty for the 
CC value is larger. It is worth recalling that this favor- 
able property of the Brandt model is counterbalanced by 
its limitation to regimes with not too high creep rates, 
whereas the CC approach is able to capture both the low- 
and the high-creep regimes. 
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